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^S^ The natural supersymmetry (SUSY) requires that stop, sbottom, and gluino be 

Oh 
-L around one TeV or hghter. By using the direct SUSY search bounds from both 

r~| ATLAS and CMS Cohaborations, we examine the constraints on the natural SUSY 

in the Next to Minimal Supersymmetric Standard Model (NMSSM). We consider 

J^ two cases of interpretations for the Higgs boson data: (1) the Standard Model (SM) 

^_i like Higgs boson is the lightest CP-even Higgs boson; (2) the SM like Higgs boson is 
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the second lightest CP-even Higgs boson. We find that the direct SUSY searches at 
the LHC impose a strong constraint on the light gluino scenarios, and in both cases 
the gluino can not be lighter than 1.1 TeV with sbottom mass up to 600 GeV and 
stop mass up to 550 GeV. 

PACS numbers: 12.60.Jv,14.80.Da,14.80.Ly 



I. INTRODUCTION 

As a leading candidate for new physics at the TeV scale, supersymmetry (SUSY) is 
strongly motivated by solving the quadratic divergence of the Standard Model (SM) and 
the gauge hierarchy problem as well as by providing a dark matter candidate and a radiative 
electroweak symmetry breaking mechanism. Compared with the situation before the start- 
up of the LHC, the discovery of a Higgs boson [H |2] and the significant constraints from the 
direct search of the LHC have driven a drastic paradigm shift in the landscape of low energy 
supersymmetry ^ [^ . We now know that the first two generation squarks must be heavier 
than 1.5 TeV or higher, though it is too early to claim the death of low energy SUSY. It 
should be noted that the models with light third generation squarks and/or light gluino, 
such as the natural SUSY models p], non-universal gaugino models [6J, and compressed 
SUSY models [7], can still be consistent with experimental data. 

The Higgs boson discovered by both the ATLAS [Ij and CMS [2J Collaborations can 
impose significant constraint on some SUSY models. For example, the naive GMSB and 
AMSB models may not produce a SM-like Higgs boson with mass around 125 GeV [S] unless 
stops are very heavy due to the small trilinear soft At term (For solutions, see Refs. P 
[To].). Meanwhile, the constrained Minimal Supersymmetric Standard Model (MSSM) or 
the minimal supergravity (mSUGRA) model may be plagued with the fine-tuning issue in 
order to accommodate the Higgs boson mass [llj via loop- induced contributions. The so- 
called fine-tuning issue can be greatly alleviated in the next-to MSSM (NMSSM) |l2j by 
utilizing both the tree-level free parameters and loop contributions [T314T7] . 

In addition to alleviate the fine-tuning issue in the MSSM, the NMSSM is also well- 
motivated by solving the /i-problem in MSSM (for a review see [E]). The singlino can ease 
the tension of experimental data and SUSY models (see [19] for more detailed discussion) 
in dark matter search. The NMSSM can be embedded into more fundamental theories, say 
the F-theory Grand Unified Theories (GUTs). The F-theory GUTs can induced unified 
boundary conditions to free parameters of the NMSSM and yield interesting low energy 
phenomenologies, as explored in [20] . 

The existence of an extra SM singlet in the NMSSM can lead to a richer Higgs phe- 
nomenology when compared with that of the MSSM, there are three CP-even Higgs bosons, 
and two CP-odd Higgs bosons. It is noticed that both two light CP-even neutral Higgs 



bosons (labelled as Hi and H2., respectively) can be SM-like. As pointed out in [21], the 
NMSSM can accommodate Higgs boson data quite well if the Higgs boson data observed 
by the ATLAS and CMS Collaborations were contributed by two degenerate Higgs bosons 
which have mass around 126 GeV. Moreover, the NMSSM has a better chance to interpret 
an extra 98 GeV Higgs boson or an extra 136 GeV Higgs boson hinted by the LEP data 
[22] or the Tevatron data [23j. It is well-known that the parameter space of either Hi or 
H2 being SM-like is different from each other ^, which motivates us to investigate these two 
cases separately. 

The null results of SUSY search at the LHC significantly and meaningfully constrain the 
masses of colored particles of SUSY [26l [27]. The experimental groups usually present their 
results in the CMSSM and some simplified models without considering any other physical 
constraints. Recently, quite a few efforts have been devoted to interpret and the LHC search 
bounds on both the MSSM [28ti3T] and the NMSSM [HI |32] with full low energy physical 
constraints. In these works, the first two generation squarks have usually been excluded up 
to about 1.4 ~ 1.5 TeV. Nevertheless, the bounds on the third generation squarks are weaker 
due to their small production rates [191 EH]. The gluino, which is have already excluded up 
to 1.4 TeV in some constrained model, can be as light as ~ 500 GeV, if the mass spectrum 
is compressed [Ml [35]. All this features are consistent with the natural SUSY spectrum, 
which motivate us to examine various scenarios of the natural SUSY, where the the third 
generation squarks and gluino may be light. 

The light gluino can play important roles in the radiative electroweak symmetry breaking, 
dark matter relic density, and gauge coupling unification at high scale. Light gluino scenarios 
can be well-motivated by Grand Unified Theories (GUTs) and string models, for example, 
the intersecting D-brane models [36H38], the F-theory GUTs [391 [IQ], the G2-MSSM [H], 
the unnatural SUSY [121 gS], the Split SUSY scenario [15H19] as well as the PeV (Split) 
SUSY scenario [SO], the natural SUSY proposed in [S], the Hidden SUSY scenario [H], and 
the compressed SUSY scenario [7]. 

Due to its large production rate, the light gluino scenarios have been a focus of phe- 
nomenological researches. Its discovery potential at the early LHC runs has been explored 



^ While in the NMSSM, the case that the heaviest CP event Higgs (i/3) might be the discovered Higgs 
boson has been explored in Rcfs. |24l[25]. 



in literatures. For example, in the reference [51], the signature of a long lived gluino under 
the split SUSY has been explored. As shown in reference [S2] , a broad and diverse sample of 
light gluino scenarios (from 350 GeV to 700 GeV) in minimal and non-minimal supergravity 
models are proposed. In the reference [53] , the scenarios of nearly degenerate gaugino masses 
are considered. The pair production of light gluino can have multi-top final states [5^ and 
multi-b final states jES], and it is expected that the multi-lepton and multi-b jet channels 
are sensitive to light gluino mass region due to the clean SM background. 

The lightest supersymmetric particle (LSP) like neutralino can be a cold dark matter 
candidate. The light gluino scenarios can also address the dark matter relic density in our 
Universe via the gluino-LSP coannihilation, which motivate us to examine such scenarios 
with current LHC data. Among them as given in Ref. [521 ES], the gluino-LSP coannihila- 
tion region is represented by points LG3-5. Such a scenario might lead to a long lifetime 
gluino and chargino, which can have interesting LHC phenomenologies, like the displaced 
kink appeared in the detectors, as explored in reference ^7] where the signature of gluino 
decaying to wino-like LSP is considered. The long lived charged wino can also be captured 
by the silicon tracker detectors at the ATLAS and CMS experiments. More recent work on 
the detection of gluino-LSP coannihilaiton can be found in Ref. |35j , where the search for 
compressed SUSY by using monojet signature has been carefully evaluated. 

Light gluino scenarios accompanied with a light third generation of squarks have been the 
intense search focus at the LHC. According to the SUSY search results from the ATLAS [2Z] 
and CMS [2S] Collaborations, gluino mass has been excluded up to ~ 1.3 TeV with rriq ~ rrig 
in the CMSSM/mSUGRA. The upper limit reduces to ~ 750 GeV with decoupled squarks. 
If interpreted in simplified models, different decay patterns have been considered separately. 
For the 100% decay chain: g — )■ tixl (mediated by a virtual i), the allowed gluino mass 
has been pushed up to ~ 1.2 TeV for the LSP lighter than ~ 400 GeV, and ~ 750 GeV 
for all available LSP masses. As for the decay chain g — t- bbxi (mediated by a virtual b), 
rrig < 1.2 TeV has been ruled out with m^o < 500 GeV. A gluino lighter than ~ 1 TeV 
with m^o < 400 GeV are excluded in the non-b tagging analysis. Bounds for a long-lived g 
which can form a R-hadron are also available, where rrig < 1 TeV are excluded by using the 
signature of slow-moving objects (low /3, fS'y) in the detectors. 

It should be noticed that all experimental bounds at the LHC are obtained by using 
some simple assumptions, where typically the decay branching fraction is oversimplified to 



be either vanishing or 100%. In reahty, to evaluate whether these hght gluino models are still 
alive or dead, model dependence must be carefully examined. For such a purpose, a more 
reliable approach is the Monte Carlo simulation, where model dependence can be correctly 
accounted for. 

In this work, we focus on the bounds to light gluino scenarios in a concrete model — the 
NMSSM. To incorporate the direct search for SUSY, we assume that squarks of the first two 
generations are heavier than 1.5 TeV, while allow the squarks of the third generation and 
gluino to be light. We explore two cases of interpretations for the Higgs boson data: (1) 
the lightest CP-even Higgs boson is around 125 ~ 127 GeV; (2) the second lightest CP-even 
Higgs boson is around 125 ~ 127 GeV. 

We find that the gluino-LSP coannihilation region in the natural SUSY models has been 
ruled out, and the gluino mass must be larger than 1.1 TeV. The reason can be attributed 
to the fact that the coannihilation region demands that the gluino is around 300 GeV. 
But, due to its huge cross section, the bounds derived from the associated mono-jet process 
pp — )■ j + gg and the Ut analysis approach require that the gluino mass be larger than 
400 GeV. It is interesting to note that the Ut analysis approach is also sensitive to the 
coannihilation region due to a relatively large fraction of boosted data sample in the signal 
events. Therefore, in the natural NMSSM, in order to be consistent with all experimental 
data, the current bounds on gluino mass must be larger than 1100 GeV (with sbottom mass 
up to 600 GeV and stop mass up to 550 GeV). 

The paper is organized as follows. In Section II, we describe the Markov Chain Monte 
Carlo method and our scanning strategy as well as the experimental bounds except the 
direct SUSY search bounds. In Section HI, we tabulate the main direct SUSY search bounds 
considered in this work, describe our workflow, and present our main numerical analysis. 
In Section IV, we examine the constraints on the benchmark points proposed in literatures 
and propose a few new benchmark points for future LHC runs and future colliders. We end 
this paper with discussions and conclusion. 



II. SCANNING STRATEGY 

A. The setup for scanning and the MCMC method 

To be general, we scan the parameter space of the natural NMSSM defined at the elec- 
troweak scale. The null search results of the signature of colored squarks at the LHC con- 
strain the first two generation squark masses to be heavier than 1.5 TeV, which motivates 
us to set their mass parameters as follows: 

MqI, = M^-^ = M^-^ = 1.5 TeV (1) 

Then we are left with a 15 dimensional parameter space (For the sake of simplicity, we fix 
Ae = 0.) to be considered . To capture the typical features of natural SUSY and light 
gluino scenarios, we choose the range of these parameters as: 

0< A < 0.7,0 < «;< 0.7,1.1 < tan /3 < 30, 100 GeV < /x < 800 GeV, 

\Ax\<3 TeV, \A^\< 500 GeV, 100 GeV < mQ^,mu^ < 700 GeV, 

100 GeV < ruDs < 1000 GeV, |A| < 5 TeV, \Ab\ < 3 TeV, 

10 GeV < Ml < 1 TeV, 100 GeV < Ma < 1 TeV, 200 GeV < M3 < 1.3 TeV 

100 GeV <Mi = M^< 500 GeV. (2) 

Obviously, it is very difficult to find the typical features of such parameter space since 
its dimensionality is too high to scan in a grid method or in a random scanning method. 
Supposed we need 10 points in one dimensional parameter space, then we need 10^^ points at 
least for the current situation. To circumvent such an issue of computational cost, instead, 
we adopt the Markov Chain Monte Carlo (MCMC) method ^5] in our scanning. 

The MCMC method is a sampling method to generate a chain of points of the parameter 
space with a density distribution consistent with experimental constraints. In this method, 
the computational time at its best is proportional to the dimensionality in a linear way, in 
contrast to the power way in the random scanning method. The MCMC method has been 
widely adopted in numerical analysis of many research fields, especially in astrophysics jSH] ■ 
Recently, this method has also been applied in SUSY model scanning [5UH53] . 

The method is inspired from the Bayesian theorem which can be put as below: 



where, p{H\d) is defined as the posterior probability of the hypothesis after taking into 
account experimental data. p{d\H) is the sampling distribution of the data assuming the 
hypothesis is true. And if it is considered as a function of the hypothesis for fixed data, it 
is called the likelihood function. p{H) is the prior probability which represents our state of 
knowledge before seeing the data. p{d) is the marginal likelihood which is just a normalization 
factor in our case and will be ignored from now. 

In this work, d and H denote a set of computed experimental observables in the model and 
a set of model parameters given in Eq. ([2]), respectively. p(H\d) is the desired distribution in 
the parameter space of our scanning after taking into account all experimental constraints. 
p{H) is taken as fiat distribution in this work, which is defined as: 

I -L (^min \ U ^ ('"max , . 

p{H) = I (4) 

I otherwise . 

p{d\H) denotes the likelihood function determined by experimental constraints, which is 
defined as 

pid\H) = l[ pid,\H), (5) 

i 

where p{di\H) denotes the likelihood function of each of experimental constraints. 

Main experimental constraints considered in this work are listed in Table [!} When the 
XenonlOO results are applied in our scanning, the proton-DM scattering cross section is 
rescaled by the formula a'p^ x Qh'^/0.11. And the likelihood functions p{di\H) adopted in 
this work can be classified into three categories: 

• Exclusion bounds imposed by setting the likelihood to be zero if the point is already 
excluded and to be one if it is not excluded. For example, to impose a likelihood value 
for all theoretical points in our scanning, we assign a zero value to a point if it is 
unphysical or theoretically unacceptable. Here a point is unphysical or theoretically 
unacceptable have quite a few meanings: the point might run into a Landau pole for 
some of its physics parameters at some energy scales, or might lead to a unphysical 
global minimum at the electroweak symmetry scale, or might have a spectra with a 
tachyonic mass for a particle, or might have a spectra of which the lightest neutralino 
is not LSP, or might fail to reach a convergent RGE solution, or might possess no 
electroweak symmetry breaking. And the bound of the XenonlOO and the bounds 
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from Tevatron and LEP on the masses of sparticles and Higgs boson are realized in 
this way. 



Upper bounds described by step functions: 

p{di\H) = 



1 



exp 



■di[H]-d^ppc 



7y (6) 

0.01 duppcr J 

where di[H] means the observable computed in the NMSSM. For example, the bound 
of dark matter relics density is realized by a step function. Although we also use a 
step function to the rare decay Bs — ?■ /i^yU^, we notice that the LHCb collaboration 
claimed a discovery of this mode with measured branching fraction 3.2]'^j^'2 x 10^^ [6^ . 
Nonetheless, our main results are not sensitive to this bound. 

Physics constraints described by Gaussian functions with well measured central values 
and deviations 



p{di\H) = exp[ 



{d,[H] 



d V 

"ccn J 1 



a^ 



(7) 



For example, the likelihood function of Higgs boson mass is taken as a Gaussian 
function with a central value 125 GeV and an allowed deviation is taken as 2 GeV. 



Experimental obserbables 


Mean vaule 


Deviation 


Ref. 


BR{B+ -^ T+Ur) 

BR{B -^ Xs-f) 


1.67 X 10-^ 
3.52 X 10-^ 


0.4 X 10-^ 
0.3 X 10-^ 


[65J 
166J 


BR{Bs -^ /u+/i-) < 4.5 X 10-9 


|67| 


f7/i2 < 0.136 


[68] 


XenonlOO (2012) 


m 


'muiggs 

Rvv 


125 GeV 
1.6 
1.0 


2 GeV 
0.4 
0.2 


imin] 



TABLE I: Physical bounds which have been taken into account in our scanning are listed here. 



Once the likelihood function, which has incorporated all experimental constraints ap- 
propriately, has been specified, we can construct Markov chains through the Metropolis 
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algorithm. The chain is a set of points in our parameter space, which can be labelled 
as {Pq, Pi, P2, ■ ■ ■ Pi, Pi+i ■ ■ ■}■ We start with a seed (labelled as Pq). The chain can be 
generated by the following steps: 1) For any a point Pi, i = 0,1, 2,---, we compute its 
value of likelihood function P{d\Pi). 2) A proposed point in our parameter space labelled 
as Pp is introduced and its value of the likelihood function is evaluated as P{d\Pp). If 
P{d\Pp) > P{d\Pi), this walk is accepted and a new start point is found Pj+i = Pp. If 
P{d\Pp) < P{d\Pi), this walk is accepted with a probability P{d\Pp)/P{d\Pi). When this 
proposed step is accepted, label it as Pj+i = Pp] when this proposed step is not accepted, 
the old point Pi will be used, i.e. Pj+i = Pj. 3) Repeat these two steps, after a long and 
sufficient walk (proportional to the dimensionality of parameter space, say 2.5 million for 
each case), a Markov Chain with sample points reflected the experimental constraints can 
be constructed. 

As observed in [60], to reduce the computational time while maintaining the sufficient 
sample points which capture the features of the constrained parameter space can be balanced 
by utilizing an appropriate proposal step in each a walk. According to the rule of thumb, 
a step with an acceptance rate around 25% ~ 30% is the best one, which is realized in our 
scanning by trial and error. 

B. Features of the Sampled Points 

In our scanning, we consider two cases of the interpretations for the Higgs boson data: 
in the first case, we assume that the Hi is the SM-like Higgs boson; in the second case, 
we assume that the H2 is the SM-like Higgs boson. We implement the MCMC method in 
NMSSMtools 3.2.1 [721 E3] and construct a chain with a 2.5 million points for each case. 
The distribution of mass spectra of sparticles for both the 1st case and the 2nd case is shown 
in Fig. [1] 

There are a few comments on Fig. ([I]) in order: 

• For the Higgs sector, in the 1st case, the mass of if 2 can spread in a quite large range 
from 120 GeV up to 1200 GeV. Similarly H3, A2 and H^ are quite heavy and also 
expand in a large range from 300 GeV up to 1600 GeV. In contrast, in the 2nd case, 
the mass of Hi is confined to be smaller than H2, therefore quite a fraction of if 3, A2 
and H^ situate below 1000 GeV. 
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Sparticle Masses for H is SM-like 



Sparticle Masses for H is SM-like 
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FIG. 1: The distribution of mass spectra of sparticles are shown for the first and second cases. 



• For the neutralino sector, in the first case, the hghtest neutrahno can spread from a 
few GeV to 340 GeV, and most of them are smaller than 200 GeV. The second and 
third neutralinos can expand from a few ten GeV to 700 GeV. In contrast, in the 2nd 
case, the LSP is compressed in a much smaller mass range from a few GeV to 220 
GeV, while most of them situate near 100 GeV. The second and third neutralinos are 
also compressed in much smaller ranges. 

• For the chargino sector, in the first case, the lighter chargino can expand from 100 
GeV to 500 GeV. In contrast, the lighter chargino can only expand from 100 GeV to 
300 GeV in the second case. 

• It is interesting to noticed that the distribution of stop squarks, sbottom squarks and 
gluino is similar in both cases. It is remarkable that in the second case, the gluino-LSP 
coannihilation region is not allowed, while for the first case, such region is possible. 

We also show the features of these points when projected on the mi^ — m^o plane, the 
mi^ — m^o plane, and the nig — m^o plane, as shown in Figs. (|2||4|). One general feature 
shown in Figs. pH) is that the favored mass of LSP is around 100 GeV for both cases. In 
the first case, the LSP is Bino or Wino dominant and its mass is determined by parameters 
Ml and M2, which is similar to the case of the MSSM. In the second case, the LSP is either 
higgsino or singlino dominant and its mass range is determined by /i. A, k, etc. 

From Fig. te|, we can read that both cases favor a relative heavy ii (say around 750 
GeV in the first case and 700 GeV in the second case), such a tendency is determined by 
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FIG. 2: The distribution of points in the m^ — m^o plane are shown for the first and second cases. 



the Higgs boson mass. 

It is remarkable that the second case allows a narrower LSP mass range than the first 
case. Such a feature can be understood by a correlation existed between the singlet scalar 
mass and the singlino mass. The mass formulae of them are provided below: 



Mhs = A\Ma ^^!! + 4k^s^ + A^Ks 



2/i 



Mg = 2ks 



(8) 
(9) 



In the second case, it is required that the lightest Higgs boson be mainly singlet-like, 
which will set an upper limit for Mh^s ^ 125 GeV. To guarantee the second Higgs boson 
mass to be 171^2 ~ 125 GeV, a cancellation condition given below 



1 - {Ax/2fL + k/X) sin 2/3 ^ , 



(10) 



is demanded, as pointed out in [13]. With this condition, from Eq. (|8]) we obtain an 
inequality 



Ak\s^ + A^Ks + X^v\l - 



K sin 2/3, 
A ' 



< 125^ . 



From this inequality we can solve out an Smax which is given as 

1 



tvSr 



A^\ + 'y/5002 + AI- 16A2i;2(l - ^1^^)) 



(11) 



(12) 



since (1 — ^'^'^ ^ ) > and most of \Ai^\ can not be larger than 300 GeV in our scanning 
(after imposing all experimental cuts), ns^ax should be smaller than 110 GeV. The only 
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exception occurs when Mh,s ^ 125 GeV due to the large mixing between the singlet and 
the doublet Hd- Then the singlet can have a mass M^^s > 125 GeV and become less than 
125 GeV after a large mixing with if^. However, this kind of space needs a certain degree 
of fine-tuning. 

Moreover, such a feature can also be understood in that a relatively small /ig// (~ 100 — 
200GeV) are favored to produce an appropriate S-H mixing so as to yield the right SM-like 
Higgs boson mass. Consequently, we have Mg = 2ns = 2^/i < /x, which result in a light LSP 
within a narrower range. 

In contrast, for the first case, the LSP is Bino or Wino dominant and its mass range is 
simply determined by the range of parameters Mi and M2 and there is no such a correlation 
exists. It is noticed that in the first case, the stop-LSP coannihilation region is allowed, 
while in the second case, such a region is missing. 




100 200 300 400 500 600 700 800 900 1000 

Mg (GeV) 



200 300 400 500 600 700 800 900 1000 

M^ (GeV) 



FIG. 3: The distribution of points in the m^ — tti^o plane are shown for the first and second cases. 



In Fig. pj), it is obvious that the first case allows more points for the sbottom-neutralino 
coannihilation regions. But generally speaking, the Higgs boson mass does not affect the 
distribution of mass of 61, since the contribution of 61 to the Higgs boson mass is much 
smaller when compared with that of ti when tan/3 < 30. 

In Fig. (li]), it is noticed that since the gluino can contribute to the Higgs boson mass via 
two loops, the mass of gluino can be affected by the Higgs boson data. The most favored 
gluino mass is around 1.1 TeV and 1.2 TeV for the two cases, respectively. 

To appreciate how stringent the LHC bounds can reach, in Figs. pH, we deliberately 
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FIG. 4: The distribution of points in the tti^^ """^v" plane are shown for the first and second cases. 



show the combined LHC bounds derived from the simphfied models. In these Figs., we 
combined the available bounds from both ATLAS and CMS collaborations |145j . In Fig. 
te|, we use the bounds of pp — > it with i — )■ txi from the ATLAS analysis based on the 
dataset of 4.7 fb~^ and ^/s = 7 TeV. In Fig. (|3J), we adopt the bounds of pp — )■ bb with 
b — !■ bxl from the CMS analysis based on the dataset of 4.98 fb~^ and ^/s = 7 TeV by using 
the aT variable. In Fig. Q, we compile two types of bounds: 1) pp -^ gg with g —^ tixi 
form the ATLAS collaborations based on the dataset of 5.8 fb~ and ^/s = 8 TeV and 2) 
g — ;■ bbxi from the CMS collaborations based on the dataset of 4.98 fb^^ and ^/s = 7 TeV. 
One obvious concern is whether all those allowed points within the bound curves are still 
alive or whether all points outside the bound curves are safe. To address this question, we 
sample 2400 points from these 5 million points (including both two cases of interpretations 
of Higgs boson data) to perform an analysis of the constraints from the direct SUSY searches 
by the LHC experiments. We have examined that the distributions of these 2400 points are 
similar to those 5 million points, therefore the analysis of these 2400 points can be reliably 
extrapolated to the sample of 5 million points. 
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III. CONSTRAINTS FROM DIRECT SUSY SEARCHES AT THE LHC 



A. SUSY search bounds from LHC experiments 



To study the constraints from the direct SUSY searches by the LHC collaborations, we 
implement in our analysis both results of ATLAS and CMS collaborations for the datasets 
accumulated with the collision energy ^/s = 7 TeV and ^/s = 8 TeV, shown in Table (UM), 
where the direct SUSY search channels from both ATLAS and CMS collaborations and the 
references are tabulated. Below we briefly describe these search channels and their sensitivity 
to possible SUSY signals: 



channels 


ATLAS 


CMS 


jets + f:^ 


mm,m\' 


um 


multi-jets + ^ 


mm,mi' 


M 


B-Jets + jets + $j^ 


[ESHEI] 


[m,mi' 


B-Jets + leptons + jets + -^ 


[EEHSa 


[m,ms\' 


letpons + jets + ^ 




[m] 


[9HHl()6j. [1171 1118)8 


multi-letpons + ^ 




mM 


|1()7H113| 


Z-boson + jets + ^ 




m[ 


Mono-jet + ^ 


[1141 


m\ 



TABLE IL The direct SUSY search results from ATLAS and CMS collaborations are tabulated, 
where the upper script 8 in the channels denotes the results being obtained from the dataset with 
the collision energy y/s = 8 TeV. 



The jets +f^ channel is the classical search channel for the signature of pair-production 
of squarks and gluino. In this search channel, a large HT{def\.ned as the scalar sum 
of pt of all reconstructed jets in one event) and a large missing transverse energy 
(denoted as pr) are required. The hard Ht is expected if all jets are produced from 
the heavy SUSY particles decay. And a large pr predicted in R-parity conserving 
SUSY can efficiently suppress QCD background. 

Apart from these two simple cuts, some characteristic kinematic variables, such as the 
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aT variable, the Mt2 variable, and the Razor variable, are used to discriminate SUSY 
signal from the SM background. Below we describe these kinematic variables in order. 

As introduced in [129j . the kinematic variable ay is designed to distinguish the real ^ 
from the hard process and the pseudo-^ from mis-measurement of jet. It's defined 
as: 



-ph ^ ^ ^ 

^ \ i=l i=l j=l 

for any a 2-jets final state (a multijet final state can be regrouped into a two-jet final 
state by using the combination algorithm that minimises the Et difference between 
the two pseudo-jets |148j ). where E^ denotes the fr of less energetic jet. A pr 
from SUSY particle decay favors a ay with values greater than 0.5, while a ^ from 
mis-measurement of jet energy typically leads to a a^ with values smaller than 0.5. 

The kink variable Mt2 is introduced in jl30j , it is supposed to determine the transverse 
mass of a new particle from its pair production with each particle decaying to a visible 
daughter and an invisible one. It is expected that two reconstructed transverse masses 
of each particle in each an event should be the same or close to each other. Similar 
to the ax variable, a multijet final state can be regrouped into a two-jet final state 
by using the hemisphere algorithm. Typically, SUSY signals can have a larger Mt2 
around several hundred GeV, while the background of the SM favor a smaller Mt2- 

The Razor variable is introduced in |131] . it is defined by CMS collaboration as: 
where M^ and Mji are defined as: 



Mn ^ J{E,^ + E,2Y - (pi' + pf )2 , (15) 



respectively. This variable has been used to search for SUSY signals with colored 
sparticles in pair production and decaying into invisible particle and jets. Signal 
events are characterized by a large Mr and a large R (which peaks around 0.5 while 
QCD multijet background events peaks around zero). In this study, we have not taken 
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into account the bounds obtained from the Razor approach, which will be included in 
our future work. 

The multi-jets + ]j^ channel is well-motivated by the signals of light squarks of third 
generation and signals of gluino decaying to squarks of third generation. Typically, 
such a signal can lead to more energetic jets, when compared with the signal of squarks 
of the first two generations. For example, the signal from pp ^ gg ^ (^^Xi)(^^Xi) with 
hadronic top quark decays can yield lots of jets in the final state. This search channel 
should be sensitive to such type of signal. 

The B-Jets + jets + ^ channel can utilize the b-tagging technique, which can be very 
powerful to reject QCD background. This search channel can improve the sensitivity 
to signals of third-family squarks production and signals of gluino pair-production with 
gluinos decaying to the third generation squarks (both on-shell and off-shell) which 
can prdocue many b jets in the final states. 

The B-Jets -|- leptons + jets + pp channel can utilize both the b-tagging technique 
and the lepton(s) (single one or two same-signed dileptons) and can reliably suppress 
the huge QCD background. It is supposed to be sensitive to the signal of four top 
final states from gluino pair production with g —> ti + ^ and pp -^ it ^ ti + ^. 

The leptons + jets + -^^ channel can utilize the high efficiency of lepton identifica- 
tion and significantly reject the QCD background. One single lepton, oppsite signed 
di-lepton and same singed dilepton channels have been considered by experimental 
collaborations. The channel is expected to be sensitive to pp — )■ titi and multi-top 
final states. 

For the multi-letpons + flp search channel, three or more well isolated leptons are re- 
quired. Trilepton channel would be golden channel to explore chargino and neutralino 
pair-production which decay to LSP and leptons mediated by sleptons. In our scanned 
parameter regions, we observe that the chargino and neutralino can be very light, as 
shown in Fig. [T| Consequently, their production rate can be very large and should be 
considered. 

The Z-Boson -|- Jets -|- ^ search channel is supposed to utilize two isolated leptons 
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from a Z boson decay. Both the momenta, sign and flavor of each leptons can be 
measured quite weU. The Z boson peak can be rehably reconstructed. This channel is 
designed to explore the topology with Z boson produced through neutralino decaying 
in the cascade decay of colored sparticles. 

• The mono-jet search channel focuses on one single energetic jet originated from the 
initial state radiation. Typically, the pt of jet is required to be larger than ~ 100 GeV. 
And a large missing energy (^ > 200 GeV) is required. This search channel can be 
sensitive to those co-annihilation scenarios, where the NLSP is almost degenerate with 
the LSP. 

ATLAS collaboration have provided upper limits for new physics in their documents and 
we use those upper limits directly. While similar upper limits are missing in the documents 
of CMS collaboration. To extract these upper limits from the CMS collaboration, we use 
the method proposed in |132[ I133J by assuming there is a 30% uncertainty on the possible 
new physics signal. 

B. SUSY Experimental Bounds Implement 

Below we will outline the main procedure as how to implement the SUSY experimental 
bounds in our study. 

For each a point selected from the constructed Markov chains, we use the NMSSM- 
tools3.2.1 [72l |73] to generate its mass spectra and decay table in SLHA format. The mass 
spectra are used to evaluate the cross sections of SUSY signals. For all points in our work, 
the most important processes include pp — ?■ gg, pp — )■ titi, pp — t- hihi ^, pp — )■ XiXj (where Xi 
include both neutralinos and charginos) pair production. We notice typically that the cross 
sections of pp — )• Xi'Xj; PP ~^ xtx^^ PP ^ xfxj can be significantly large due to their small 
masses, like in the second case. 

The NLO cross section is evaluated by using the package prospino2 |139j . which will be 
used to normalize the number of signal events in our analysis. Then, the mass spectra and 
decay tables are passed to the package MadGraphS |140] . and signal events with upto two 

^ We have taken into account the contribution of the process pp — ?> 6262 and have found that the LHC 

bounds start to constrain those cases with Mr < 600 GeV. 

02 
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additional radiative jets for the processes pp — > gg, pp — > titi, pp — )■ 6161, pp — )■ xtXj are 
generated. To avoid double counting issue in the matrix element calculation and the parton 
shower simulation, we adopt the MLM-matching scheme with the variable xqcut = lOOGe^. 
Then PythiaG |141j is used to decay the sparticles to the particles of the SM at parton 
level and to simulate parton shower and hadronization. We use PGS4 [142J to implement 
fast detector simulation. To reconstruct jets in the final objects, we adopt the anti-kt jet 
algorithm with the cone size parameter R = 0.5, and assume the b-tag efficiency to be 60% 
in accompany with a mis-tagged rate for charm quark jet as 10% and for other light quark 
jet as 1%, respectively. 

We generate 50,000 events for each of signal processes at parton level, after matching 
typically we arrive at 30,000 matched events or more. The matched events will be passed 
to our SUSY bound analysis package to evaluate how many number of events can survive 
after implementing all experimental cuts. 

To analyze the bounds imposed by the direct SUSY search at the LHC, we develop 
a systematic analysis package. The main goal of the package is to implement the SUSY 
constraints given by ATLAS and CMS collaborations in an ease way. With the help of our 
package, we can evaluate whether a model is still alive or has been ruled out. Up to now, 
the searches that we have implemented are given in the Table |TT| and we are upgrading our 
package by including new LHC bounds released recently. 

For each search channel by feeding the matched events of signal processes to our package, 
we can read out the selection efficiency for in each signal region. This selection efficiency 
finally is translated into the observed number of signal events in each signal region after 
cross sections and luminosity are taken into account. In order to perform an analysis similar 
to the LHC collaborations, for each a point we generate two independent event samples with 
the collision energy y/s = 7 TeV and y/s = 8 TeV, respectively. 

In this work, the exclusion limits up to the observed 95% confidence level for each search 
channel in each signal region have been applied. Accordingly, we define the ratio R = 
- for each signal region in each searching channel. To derive the most stringent 



-'*signal iiumbc 

-'^observed limit 



constraint, we choose the maximal value of R among all search channels at all signal regions. 
Obviously, for a model at a specific search channel and a specific signal region, the ratio R 
is greater than 1 means it has been ruled out by experiments (although we have not taken 
into account errors, neither Monte Carlo errors nor the fast detector simulation errors). Our 
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package yields all Rs of each signal region at each search channel. By comparing the Rs, we 
can find out the strongest bound, which is denoted by -Rmax- 

C. Numerical Analysis 

Here we present our main results of numerical analysis. In Fig. (|5|, we present the cross 
sections of gluino and stop in our scanning. We highlight two observations: 1) Roughly 
speaking, the cross sections of stop and gluino increase by a factor two when the collision 
energy increases from 7 TeV to 8 TeV; 2) when stop and glunio have the same mass, the 
cross section of gluino is 50 times larger than that of stop. It is also noticed that the cross 
sections are almost equal for the stop pair production and the gluino pair production if the 
mass of 171}: ^ rrig — 250 when rrig = 600 GeV and mi ^ rrig — 400 GeV when rrig = 1200 
GeV. At the tree-level, the cross sections of pp — )■ tt and pp -^ gg are simply determined 
by the mass parameters, while at the NLO level, colored sparticles at loop can contribute 
and lead to a minor change in the cross section. It is noticed that the fiuctuation in the 
cross section of the process pp — )■ tt near the mass region 500 ~ 600 GeV by a few points is 
originated by the stop-decay threshold {mi = mt + mg) effect [TMt 1135] . The cross sections 
of pp — )■ hh is similar to that of pp -^ tt. Therefore we neglect them in the Fig. (pj). 

In Fig. (pi), we first examine the constraints on the signal of either the process pp — )■ tiii 
or the process pp — ?■ hibi in the m — m^o planes, as shown in the right plot and in the left 
plot, respectively. The experimental bounds is the same as shown in Fig. ([2|. The Rmaxiii) 
or Rras,x{bb) is obtained by using all kinds of SUSY search analysis approaches (including 
both ax and Mt2, etc.). 

It is interesting to notice that the green points inside the bound curves simple indicate 
that the branching fraction is too small and yield too small number of signal events to be 
meaningfully constrained. While the black points outside the bound curves are found to 
be constrained by other search channels not deliberately designed for either the signal of 
pp — )■ iiii or the signal oi pp — )■ bibi. For example, the black points outside the bound curve 
near the point [500, 100] in the m^ — -m^o plane and those near the point [600, 100] in the 
?7ij — m^o plane are ruled out by Jets (M7-2) search channel, while several points near the 
region of point [500, 250] (with b — )■ bx^, x^ — > //, and slepton decaying to lepton and LSP) 
are constrained by the 2SSL search channel. These black points outside the experimental 
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FIG. 5: The cross sections in varying with the mass of gluino and stop are shown here. The cross 
section of sbottom is similar to that of the stop and is omitted. 



bound curves clearly demonstrate the importance and necessity of a comprehensive analysis 
for a given model. 

Compared with the results given in [T^ 1136] , we observe that when more experimental 
constraints up to 5 fb^^ with ^/s = 7 TeV and part of those from the analysis with ^/s = 8 
TeV are included, the bounds to stop and sbottom have been improved, as clearly demon- 
strated by the scattering plots in the upper row of Fig. ([t]). It is straightforward to read 
out that the stop mass can be excluded up to 550 GeV or so, while the sbottom mass can 
be excluded up to 600 GeV or so. 

In Fig. (It]), we show the bounds for the stop, sbottom and gluino. In the lower row, 
we show bounds against the gluino mass. We observe that most of points have been ruled 
out /disfavored even if they can be constrained meaningfully neither by the signature of 
pp —7- titi nor by the signal of the process pp — > bibi. When gluino is light, say less than 
800 GeV, the most stringent bound is from the signatures of pp — )■ gg due to its large cross 
section and the bounds show a universal model dependence which is indicated by the width 
of the band. The width of band when rrig < 800 GeV is around 10, i.e. the sensitivity to 
signature of a given rrig can differ by a factor of 10. In contrast, the sensitivity of signature 
of pp — >■ bibi can differ even larger. 
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FIG. 6: The constraints to the sohtary signal oipp — )• titi and pp — )• 6i6i are demonstrated, respec- 
tively. -Rmax(i^)(-Rmax(&&)) means the largest R values in all search channels while the signature 
only includes the process pp — >■ titi {pp — )■ bibi). 



Comparing the left and right plots in the lower row, we observe that the bounds have a 
strong correlation with the sbottom mass. This correlation, especially at the right corner 
with 1 TeV < rrig < 1.3 TeV where lots of points with a heavy gluino have been excluded, 
can be attributed to the fact when the cross section pp — )■ gg is much less than that of 
pp — )■ bibi, consequently the real meaningful constraint is actually from pp — t- bibi. While 
such a correlation with the stop mass is weak. The width of band near the region 1 TeV 
< rrig < 1.3 TeV become broader, since it is determined by the signature of pp — > bibi, 
instead of pp ^ gg. 

In Fig. (Is]), we show four representative constraints to the gluino signals from different 
analysis approaches which are supposed to be sensitive to the signature of the production 
process pp ^ gg ^ ttttfrp, for instance. In the upper left plot, the bound from the multijet 
analysis approach is demonstrated, which can exclude the signals of most models with gluino 
mass lighter than 600 GeV. In the upper right plot, the bound from the B-Jet plus ar 
analysis can exclude the signals below 900 GeV. Meanwhile, this analysis approach enjoys 
less model dependence than the multijet analysis approach as indicated by the width of the 
band formed by the points. 

In the lower left plot, the bound from one lepton plus B-Jet plus -^ is shown and the 
meaningful constraints can reach up to 800 GeV. At the lower left corner of this plot, there 
are some points can not be constrained due to the small branching fraction oi g ^ ti + fJr 
and the dominant branching fraction is g ^ bb-\- ^ or g ^ g + pr- While in the lower right 
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FIG. 7: In the upper row, the bounds for stop and sbottom are examined where the signals are 
assumed to be pp — )• iiii and pp — )• 6161, respectively. The y-axis is the Rma.x, and the x-axis is 
the mass of stop and sbottom, respectively. In the lower row, the bounds for all signals (including 
pp — )• titi, pp — )• bibi, and pp — )• gg, are taken into account. The x-axis is the mass of gluino, the 
y-axis the Rmax, and the color bar indicates the mass of stop and sbottom, respectively. 



plot, the bound from the same sign lepton is demonstrated and the meaningful constraint can 
reach up to 800 GeV or so. Obviously, both of these two channels rely upon the branching 
fraction of ^ — )■ tt + ^. Although the same sign lepton mode is clean and has a very tiny 
SM background, its sensitivity is similar to or worse than the search channel of one-lepton 
+ B jet + prp due to its much smaller branching fraction. 

Among these four search channels, it is worthy of remarking that the bounds from the 
ot analysis approach with b-tagging is the most stringent and the least model dependent. 

We also notice that similar to the ar approach, the Mt2 approach also enjoy a good 
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FIG. 8: Four channels for the signals from pp — )• gg are shown. The x-axis is the gluino mass, while 
the color scale denotes the mass of LSP. 



sensitivity and model independence, as shown in the left plot of Fig. (tol), where both the 
constraints from the Mt2 approach and the constraints from the mono-jet plus ^ channel 
are shown. We notice that the Mt2 approach can achieve a sensitivity better than the ax 
approach for most of points, as be revealed by the statistics shown in the right plot of Fig. 



(10). 



Meanwhile, the mono-jet search channel alone can probe the gluino mass up to 400 GeV, 
as shown in the right plot of Fig. (jol). To show the constraint of the mono-jet on the 
gluino-LSP coannihilation scenarios, we deliberately introduce extra 200 points in this plot 
as denoted by the markers of an empty diamond. 



At the left plot of Fig. (10), 50 representative points with their branching fractions 



denoted by pie charts are shown to demonstrate the effects of branching fractions of four 
main decay chains (say g -^ tiflj,, g -^ bbflj^, g — )■ gflp, and others). For each category 
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FIG. 9: The constraints from the Mt2 analysis approach and the mono-jet channel are demon- 
strated, respectively. In the right plot, extra points of the gluino-LSP coannihilation are displayed 
by the markers of an empty diamond. 



of decay modes, we sum over all on-shell or off-shell decay modes into one. For example, 
g — 7- t^*H^*^ means we count on either the on-shell or off-shell decay modes oi g -^ ti, g -^ t*i, 
g —^ ti*, g -^ t*i*, and sum over all allowed decay modes. While g -^ bb^*^ means either the 
on-shell or off-shell decay modes oi g ^ bb and g ^ bb*. 

It is noticed that the most stringent bounds of gluino with mass lower than 1000 GeV do 
depend upon the branching fractions, as sensible from the width of the band of -Rmax, which 
denotes the difference of sensitivity for a given rrig. This difference can change by a factor 
of 10 when rrig < 1000 GeV, which can be attributed to the fact that when the cross section 
of pp — >■ gg is large enough and the branching fractions of ^ — t- 66 + ^^ and g —^ ti+ ^rp can 
be large enough, the constraint the Mx2 analysis method with b-tagging can always have a 
good performance, though the branching fraction of ^ — ;■ tipr does modify the sensitivity 
to some degree. When rrig > 1000 GeV, the dominant signals might from either pp — )■ bibi 
or pp — )■ titi, then the model dependence of the bounds increases. 

It is also noticed that when gluino is around 300 GeV ~ 500 GeV and dominantly goes 
to g + $j.^ due to the large cross section of pp — > gg and the large mass splitting between 
the gluino and the LSP (say larger than 100 GeV), either Mt2 or ax analysis approach 
can have a remarkable sensitivity to these points. Due to the fact that we assume that the 
squarks of the first two generations are heavier than 1.5 TeV, therefore for most of points 
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the branching fraction of (7 — )■ qq' + $j' is typically negligible. 
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FIG. 10: At the left plot, the branching fractions dependence is examined by choosing 50 represen- 
tative points. At the right plot, the statistics on how inany inodels are excluded by which channels 
are shown. 



At the right plot of Fig. ( 10 ), among these 2400 points, the statistical information on how 



many models are excluded by which channels is provided. It is easy to read out that both 
the ut and Mj'2 analysis approaches are overwhelmingly sensitive to most of points, while 
the search channels with lepton(s) (especially the multiletpon channel with the required 
lepton number larger than ni > 3, in which case the standard model background is almost 
vanishing) start to play a role when neutralinos and charginos are light and their production 
rates are large. It is also remarkable that when gluino is light and the decay mode g — )■ bbxi 
is dominant (say around 300 ~ 500 GeV), the search channels with b-tagging can be efficient. 



In Fig. (11), we project all points in the rrig — m^o plane so as to compare with the 



experimental bounds directly, which have also been shown in Fig. (El). There are quite a 
lots of points outside the experimental bounds are excluded mainly due to the light sbottom 
in the mass spectra as marked out at Fig. ([T]). 

In order to examine those points of near the gluino-LSP coannihilation region, in the first 
case, we deliberately introduce 200 extra points which are denoted by markers of an empty 



diamond near the points [300,300] and [400, 360] in Fig. (11). It is found that these points 
can be excluded by both the search channel of mono-jet + fip and the search channel jets 

It is remarkable that for these coannihilation points, the search channel of jets plus pr 
(say, B-Jets + ^, Mt2 and ax analysis approaches) has a better sensitivity than the channel 
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FIG. 11: The scatter points projected in rrig — m.yO plane are provided, where the black points are 
heavily disfavored and the green points are safe. The gluino-LSP coannihilation region of points 
[300,300] and [400, 360] are denoted by markers of an empty diamond. 

of mono-jet plus ^ . The underlying reason is that due to the large fraction of boost data 
sample in the signal events (about 3 ~ 6% of the total cross section of pp — )■ gg) when 
compared with the cross section of the mono-jet events after using the mono-jet search cuts 
(about 1 ~ 2% of the total cross section of pp — > gg). 

For the second case, there is no coannihilation region with mass splitting smaller than 20 
GeV. We notice the bound of gluino mass is close to the simplified model due to the large 
mass splitting Am = rrig — m^o and the energetic final states. 

IV. BENCHMARK POINTS 



We examine benchmark points recommended in literatures. We consider 6 benchmark 



points labelled by "NMP" in Table |n!| of NMSSM listed in [137] and find all of them have 
been excluded by experiments to a quite high confidence level. We also check two benchmark 
points labelled as "EHP" of NMSSM from pjSSJ and two benchmark points labelled as 
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"DET" of NMSSM from [32]. We also examine the so-called light slepton benchmark point 
{5Mf) compiled in Table 2 of [3] with a light slepton sector inspired by anomalous magnet 
momentum of muon of pMSSM (where both stop and sbottom are also light). For each of 
benchmark points, we perform the same analysis as those points given in Section III. The 



bounds obtained from each channel are listed in Table (III) 



We present our results in form: -Rgg^rch 



channel, signal region' 



where R is the maximum ratio 



of N^ig/N^i among all the search channels in each category of signal regions. And the letter 
"A" means those from ATLAS collaboration while the letter "C" means those from CMS 
collaborations. 
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TABLE III: The maximal ratio of r = N/Nexp in each category of all channels are shown here, 
where A^ denotes the number of events after all cuts and N^^p denotes the allowed number of events 
by experiments. Numbers in brackets are CLg for the signal. Super- and sub-scripts denote the 
LHC collaboration (A means ATLAS and C means CMS), the collision energy {^/s = 7 TeV and 
-y/s = 8 TeV) and the search channel, respectively. The label "tmp" appeared in the 5th column 
denotes the results obtained by the template approach used by the CMS collaboration |146j . 



For the original benchmark points labelled as "NMP" and proposed in |137j . we observed 
that the search for the first two generation squarks already rule out all of them. So we 
modify the mass of the first two generation squarks to 1.5 TeV. It is noticed that the first 
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six benchmark points are very similar in spectra and decays. Each of these six benchmark 
points has a gluino mass around 700~800 GeV with electroweakinos hghter than ~500 
GeV. It is noticed that for all these six benchmark points, gluino dominantly decays to tit. 
Therefore there is no surprise that the search channels, like B-Jet(s) + ^, lepton + B-Jet 
+ jets + ^, and same sign lepton + B-Jet + jets + fij,, are sensitive to the signals of these 
points. Due to the large mass splitting Am = rrig — m^o, the gluino-LSP coannihilation can 
not occur, which consequently yields energetic visible final states when gluino goes to LSP 
and consequently leads stringent constraints for all these benchmark points. 

As observed in [12], these benchmarks points can survive the light stop search bounds 
when only part of direct search bounds are applied. When more direct search bounds are 
applied, to save these six benchmark points, it is possible by assuming that gluino mass is 
higher than 1.5 TeV, as done in [19]. Except glunio mass, for the third benchmark point, 
the stop mass must be higher (say 500 GeV, for instance). 

For two benchmark points labelled as "EHP" and proposed from |138j . it is observed that 
the most stringent bound is obtained from the search channels jets+^, B-Jets + -^ (with 
ttT and Mt2 analysis approaches), and lepton + jets + flj., which are sensitive to the signal 
from the production process pp — )■ titi. These two benchmark points have been studied in 
[19] and have been found to be marginally safe when only the ATLAS analysis with 2B-Jets 
+ -^ by using 2.05 fb~^ dataset and the CMS analysis with 2B-Jets + ^ by using 4.98 fb~^ 
dataset with ^/s = 7 TeV are applied. Now with more results from the ATLAS and CMS 
analysis, these two benchmark points start to be in trouble. These two benchmark points 
can only become safe when the top squark mass is lift up to higher than 600 GeV. 

About the two benchmark points labelled as "DET" and proposed recently in [32] are still 
marginally safe. The dominant signals of these two benchmark points are from the process 
pp — )• tt, though the signals from electroweakino's pair production are also considerable. It 
is observed that the Mt2 analysis method can put even more stringent bounds that other 
analysis method. When the updated SUSY bounds are taken into account, only the first 
point is marginally safe. 

About the light slepton benchmark points labelled as "5Mf " and compiled in ^ , although 
the gluino is heavy, this point has been excluded by the LHC experiments. The most sensitive 
channel is from the multilepton search mode due to the large cross sections of slepton and 
electroweakinos' pair production, other search channels sensitive to colored objects also 
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disfavor this point due to its light stop and sbottom. 



Based on our analysis, we propose four benchmark points tabulated in Table (IV) which 
are safe and can be examined for the future LHC runs. 



Points 


I 


II 


III 


IV 


A 


0.648 


0.673 


0.349 


0.499 


K 


0.323 


0.252 


0.415 


0.140 


tan/3 


2.71 


2.68 


22 


9.8 


/^cff 


303 


509 


116 


213 


Ax 


641 


1208 


2626 


2261 


A^ 


-362 


-231 


-446 


-110 


^Qa 
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938 
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811 


^C/3 
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620 


980 


964 


^D, 


275 


763 


820 


957 


At 


1792 


1450 


1745 


1833 


Ab 


-60.6 


-2903 


2887 


2910 


Ml 


854 


530 


786 


522 


M2 


964 


269 


493 


260 


M3 


1094 


1013 


1155 


1174 


Ml 


380 


371 


212 


322 


niH^ 


125.7 


125.8 


118 


95 


mH2 


191 


335 


126.2 


126.3 


niH; 


827 


1431 


2625 


2186 


rriA^ 


424 


392 


421 


149 


rfiA^ 


824 


1429 


2625 


2186 


mu± 


819 


1425 


2625 


2184 


m-g 


1209 


1130 


1267 


1285 


m^o 


260 


260 


99 


98 


"^x^ 


-324 


381 


-129 


208 
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I 


II 


III 


IV 


m^o 
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496 


295 


-235 


^xl 


841 


-523 


518 


317 


"^xi 
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570 


775 


517 


m-± 

Xi 


299 


263 
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181 


"^xt 


979 


532 


519 


317 


^h 


746 


589 


790 


741 


^h 


1064 


993 


1081 


1063 


^b. 


332 


795 


848 


834 
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958 


878 


991 


^UL 


376 


367 


202 


316 


"^e^/ML 


382 


373 


217 
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"^efl/Afl 
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373 


216 


325 


ruf^ 


380 


370 


206 


319 


rrir^ 


383 


376 


227 


331 


BR{B+ -^ T+Ur) X 10^ 


1.32 


1.32 


1.31 


1.32 


BR{B -^ Xs-i) X 10"^ 


3.63 


3.40 


3.89 


3.56 


BR{Bs ^ n+fi-) X 10^ 


3.68 


3.68 


3.69 


3.68 


nh^ 


0.006 


0.003 


0.01 


0.10 


a^^ (pb) X 10^ 


26.6 


31.2 


3.2 


1.1 


77 


1.06 


1.02 


1.25 


1.14 




1.03 


1.0 


0.98 


1.05 


BR{g^iit){%) 


34 


65 


46 


53 


BR{g^bib){%) 


59 


27 


28 


30 


BR{g^b2b)i%) 


7 


8 


26 


15 



TABLE IV: Benchmark points for future LHC runs are tabulated. 
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For each case, two benchmark points are presented (mass of LHC Higgs is shown in 
boldface). For the first case, benchmark point I has a very hght sbottom ~ 300GeV, 
which survives from the constraints of the SUSY direct search owing to a heavy LSP(260 
GeV). Such a benchmark point might be probed by the full dataset collected with ^/s = 8 
TeV, as shown in |147] . Benchmark point II has a relatively light stop (589 GeV), which 
dominantly decay to X2 t (29%) and xt b (61%) with X2 ^ Xi W^{90%) and xt ^ 
xt Z (32%), xf Hi (23%), Xi W^ (26%), xl W^ (13%). Longer cascade decay chains 
result in more objects in the final state, many of which are too soft to be reconstructed by 
detectors of the LHC. 

For the second case, benchmark point III and IV have relatively heavier stops and sbot- 
toms which are safe to current LHC constraints. However, point IV have all five neutralinos 
and two charginos very light (wino masses are only 317 GeV). It would be expected that 
the neutralino-chargino search(mainly through tripleton and same-signed leptons signal) at 
LHC should be sensitive to this point. And this point gets correct relic density by virtue of 
the large singlino component (~ 64%) of LSP, while the other three points has LSP mostly 
Higgsino-like (I and HI) or Wino-like (II). 

It is noticed that masses of the gluinos are around ~ 1100— 1300GeV. The main branching 
fractions of dominant decay modes of gluino are tabluated. As one can see, the gluino 
will dominantly decay to tit and 61/2& which are exactly the representive simplified models 
explored by ATLAS and CMS. These benchmark points can be detected in the future LHC 
runs. And the recent SUSY search results with ^/s = 8 TeV presented in the Moroind 2013 
EW could be sensitive to these points. 

Another interesting fact is that all the sleptons are lighter than ~ 500 GeV, currently 
LHC results can not impose meaningful constraint to slepton sectors (The LHC results can 
probe the sletpons with masses up to 200 GeV). The future LHC runs at higher collision 
energies, either ^/s = 13 or ^/s = 14 TeV can start to probe them. 



V. DISCUSSIONS AND CONCLUSIONS 

We have extended the analysis [121 I136J by including more LHC experimental bounds. In 
order to test the reliability of the results, we perform a thorough check to our package. We 
have used 21 test points in total to compare with experimental results. These 21 test points 
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have been tabulated in Table (fVl), where TPl-7 denote benchmark points of the CMSSM 



as shown in Fig. 12 And GBs and GTs are benchmark points of the simplified model, 
with g — )■ tixi or ^ — )■ bbxi with a branching fraction 100%. In addition, the test points 
labelled as "EW1-EW7" denote specific benchmark points designed for the electroweakinos 
search |149j . whose decay modes are also displayed. For test points "EW1-EW5", we assume 



mr 



0.5m-± + 0.5m^o. 



Xi 



Xi 





Mo 


Ml/2 






^^^xVxt 


M-o 






M-g 


M-o 

Xl 




TPl 
TP2 
TP3 
TP4 
TP5 
TP6 
TP7 


210 
230 
85 
500 
700 
1450 
2000 


285 
360 
400 
500 
600 
175 
300 


CMSSM, LMl 
CMSSM, LM5 
CMSSM, LM6 

CMSSM, LM9 


EWl 
EW2 
EW3 
EW4 
EW5 
EW6 
EW7 


400 
500 
400 
400 
350 
200 
200 


200 
300 
250 
100 
150 
75 
50 


x'^2^ll{BF = 0.5), 

xf -^ 10, ul 


GBl 
GB2 
GB3 
GTl 
GT2 
GT3 
GT4 


800 
1000 
1050 
700 
800 
900 
1000 


300 
400 
550 
100 
150 
250 
350 


9 ^ bbx'l 


xl^ll{BF = 1), 
xt -^ i^tTr 


9 -^ tixi 


xl ^ Zxl 

Xt ^ M/±x? 



TABLE V: All test points are tabulated. 



The results of 21 points are tabulated in Table (VI). In the table, search channels are 



arranged by the order of the ATLAS and CMS Collaborations and in the order of collider 
energies, with datasets of ^/s = 7 TeV first and followed by the datasets of ^/s = 8 TeV. 
For the multilepton channel, we list the check results separately. 



In each row of the Table (VI), one can read the results for each search channel. For 



example, it is noticed that the search channel "1-2B-Jets + 1-2L" is insensitive to all our 
test points, since this channel is designed for very light stop (with mass similar to, or lighter 
than top quark) search. 



From Table (VI), we can observe a fact that for test points "TP1-TP7", "GB1-GB3" 



and "GT1-GT4", the Mt2 observable is quite sensitive to most of signals, similar to our 
observations in Section III. Furthermore, the results for the "Multilepton" search channel 
also agree with experimental ones [149 j very well. It is also obvious that the sensitivity of 
a specific search channel can perform better for the dataset of ^/s = 8 TeV than that of 
^/s = 7 TeV due to the enhancement of cross sections, as demonstrated in the channels 
"Multijet" of the ATLAS and "2SSL + B-Jets" of the CMS. 

In the "Jets+MHT" search channel, one can read out that for seven CMSSM test points 
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five are all excluded except for TP5 and TP7, which is in agreement with the results of 



CMS-7TeV as shown in Fig. 12 While the result for GBl shows a minor deviation. As 
shown in Fig. 13, this test point should have been excluded by the search channel "Jets + 
Ifrp^\ while our result underestimates the constraint with a -Rmax = 0.8. 

We observed that most of our results agree with experimental results, though there are 
some results showing deviations. For example, the "GB3" is under-constrained by the 
2SSL + B-Jets channel, and the "GT4" is over-constrained by the Mt2 analysis approach 
and is under-constrained by the multijet channel. Compared with experimental results. 



these deviations of R values in Table (VI) can be typically around ±30% or so, which 
can only affect the results of those points near the exclusion edges of experimental results. 
For example, the results of "GB3" and "GT4" should be excluded by the corresponding 
experimental search channels, but can survive in our analysis. When a point is far away 
from the edge and inside the exclusion region, our results are trustable, like the test point 
"GTl". 

Meanwhile, we notice that when the gluino mass is shifted by ±40 GeV, the deviations 
can be removed and agreements can be achieved. The deviations in R can be attributed 
to the small fluctuations from Monte Carlo simulations and the difference between the fast 
detector simulation and the real detector effects. Our bottom line is that our results at least 
motivate experimentalists to perform a more detailed and serious bounds when all types of 
uncertainties are taken into account. 

Based on our analysis shown above, below we comment on light gluino models in litera- 
tures. Models in the coannihilation region as demonstrated by LG3-5 given in [52] can be 
ruled out due to the mono- Jet search and by the ar and Mt2 analysis (due to the boosted 
sample events) from the CMS if the branching fraction g — > f^x^ is dominant. Due to the 
large mass splitting between the gluino and the LSP and the energetic final states, light 
gluino scenarios with mass range [500,700] GeV in the so-called Higgs pole region of the 
minimal universal SUGRA |143j is heavily disfavored. The two benchmark points presented 
in Ref. [55] has been excluded due to their large cross sections and large mass splittings. 
We also find that even for the coannihilation region (say nig — niyX) < 20 GeV) the mono-jet 
bounds can be valid for gluino mass up to 400 GeV at least, as shown by the right plot of 
Fig. (§. 

The last but not the least, in our scanning, we have not found points near the region 
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TP4 
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21 
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0.01 


2.7 
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TABLE VI: Results for each individual search channels are presented for checking. Numbers in 
the table are defined as R = Ngig/NuL, where R > 1 means excluded by experiments. The slots 
with a hyphen mean we haven't look at these channels since the signal is expected to vanishing. 
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rrig ^ m^o ~ 500 GeV. Such models should exist since there is no upper bounds for the LSP, 
as demonstrated in the Fig. (8c) of [144] by few red points, where the scanning ranges for 
parameters of neutralino and gluino are wider than those in our scanning. 

With results shown above, we would like to conclude: 1) the parameter space of the 
natural NMSSM has been significantly shrunken by the direct SUSY search bounds from 
the ATLAS and CMS collaborations; 2) The gluino-LSP coannihilation region in the first 
case of interpretations of the Higgs boson data has been completely ruled out by the boosted 
signal samples. We can safely claim that the lightest gluino mass must be larger than 400 
GeV for the gluino-LSP coannhilation region. 

Note: When this work is finished, in the updated results presented at the Rencontres de 
Moroind 2013 EW, we notice that the gluino mass must be heavier than L3 TeV when 
the squarks of the third generation are assumed to be light. Experimental bounds nearly 
approach the upper limit of gluino mass of the natural SUSY. 



3 production, g->ttx 



CMS Preliminary L = 4.98 fb"\\/s = 7 TeV 




500 1000 1500 2000 2500 3000 

mo [GeV] 






, 600 



200 



100 



OL3 95% C.L. limits, a^^^ not inciuded. 
0-iepton. > 6-9 jets [L = 5.8 fb"', 8 TeV] 

ATLAS-CONF-2012-103 '" 

2-SS-ieptons, > 4 jets [L = 5.8 fb"', 8 TeV] 

ATLAS-CONF-2012.105 '" 

3b-jets [L = 4.7 fb', 7 TeV] 




FIG. 12: Test points labelled as TP1-TP7 of the CMSSM in the mQ — mij2 plane and those labelled 
as GT1-GT4 in the simplified models are displayed in our comparison with experimental results. 
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FIG. 13: Test points labelled as GB1-GB3 and GT1-GT3 in the simplified models are displayed in 
our comparison with experimental results. 
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